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Patient-specific induced pluripotent stem cells (iPSCs) will 
assist research on genetic cardiac maladies if the disease 
phenotype is recapitulated in vitro. However, genetic back- 
ground variations may confound disease traits, especially 
for disorders with incomplete penetrance, such as long-QT 
syndromes (LQTS). To study the LQT2-associated C.A2987T 
(N996I) KCNH2 mutation under genetically defined condi- 
tions, we derived iPSCs from a patient carrying this muta- 
tion and corrected it. Furthermore, we introduced the same 
point mutation in human embryonic stem cells (hESCs), 
generating two genetically distinct isogenic pairs of LQTS 
and control lines. Correction of the mutation normalized 
the current (/kJ conducted by the HERG channel and the 
action potential (AP) duration in iPSC-derived cardiomyo- 
cytes (CMs) . Introduction of the same mutation reduced /xr 
and prolonged the AP duration in hESC-derived CMs. 
Further characterization of N996I-HERG pathogenesis re- 
vealed a trafficking defect. Our results demonstrated that 
the C.A2987T KCNH2 mutation is the primary cause of the 
LQTS phenotype. Precise genetic modification of pluripo- 
tent stem cells provided a physiologically and functionally 
relevant human cellular context to reveal the pathogenic 
mechanism underlying this specific disease phenotype. 
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Introduction 

Transgenic animal models of human cardiac diseases and 
cultured cells have been crucial to cardiovascular research, 
contributing to both our basic understanding of normal 
development and disease mechanisms. However, these ap- 
proaches often only partially recapitulate the molecular and 
cellular phenotypes observed in patients because human 
heart physiology, as well as the cardiomycyte's gene and 
protein expression profile, differs from almost all other 
experimentally accessible species (Davis et aU 2011). 

Human induced pluripotent stem cells (hiPSCs) (Takahashi 
et aly 2007b; Yu et al, 2007] are creating exciting new 
opportunities for biomedical research by providing 
platforms to study mechanisms of genetic disease 
pathogenesis that could lead to new therapies or reveal 
drug sensitivities (Braam et aU 2013; Liang et at, 2013). 
Patient-specific hiPSCs that carry all disease-relevant genetic 
alterations are important not only for understanding 
monogenic and complex disease mechanisms in the 
affected cell types in vitro and identifying potential 
treatments (Gold-von Simson et al, 2009; Lee et aU 2009; 
Marchetto et al, 2010; Axelrod et at, 2011; Pini et al, 2012], but 
also for providing insights into factors that predispose 
individuals to develop overt symptoms. However, the 
application of hiPSCs in disease modelling has been limited 
by the paucity of experimental tools for distinguishing subtle 
but disease-relevant phenotypic changes from background- 
related variations. 

Inherited long-QT syndrome (LQTS] is a life-threatening, 
often autosomal dominant, disorder characterized by pro- 
longed ventricular repolarization, a propensity to poly- 
morphic ventricular tachycardia and sudden cardiac death 
in young patients (Crotti et aU 2008]. Among genotyped 
individuals, ~90% of affected patients have mutations in 
the genes encoding repolarizing channels of the delayed 
rectifier currents, Jks [KCNQl in LQTl] and Jxr [KCNH2 in 
LQT2] (Morita et al, 2008]. LQTS presents clinically with a 
broad range of phenotypes even among family members with 
identical mutations (Giudicessi and Ackerman, 2013]. This 
heterogeneity is thought to result from the interplay of 
complex networks of direct and indirect causal factors, 
including epistatic effects of genetic background, resulting 
in incomplete penetrance and variable severity. LQTS was 
one of the first cardiac diseases recapitulated using hiPSCs as 
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a model (Moretti et al, 2010; Itzhaki et al 2011; Matsa et al 
2011; Yazawa et al 2011; Davis et al 2012; Egashira et al 
2012; Lahti et al 2012). However, there is a growing 
realization that comparing patient cells to unrelated healthy 
control cells, no matter how they are matched, does not 
control for any effects of genetic modifier loci (Bellin et al 
2012; Cherry and Daley, 2013; Sinnecker et al 2013). 
Furthermore, additional key factors may contribute to 
variability among individual hiPSC clones, including the 
number and location of viral integrations following 
reprogramming (Gonzalez et al 2011), genetic heterogeneity 
in the starting somatic cell population (Gore et al 2011), copy 
number variations (Panopoulos et al 2011; Abyzov et al 
2012), and their epigenetic status (Nazor et al 2012). For 
these reasons, it is preferable to compare cell lines in which 
the disease-causing genetic lesion of interest is the only 
modified variable. Gene correction in hiPSCs from patients 
with known genetic mutations is a powerful tool for 
overcoming the limitations imposed by individual 
variability among independent controls, and with recent 
technological improvements genome editing is now feasible 
in both human embryonic stem cells (hESCs) and hiPSCs 
(Costa et al 2007; Zou et al 2009; Hockemeyer et al 2011; 
Asuri et al 2012). The generation of genetically matched 
controls that differ exclusively in well-validated susceptibility 
variants for maladies is a generally appHcable solution to 
the key problem of distinguishing pathologically relevant 
phenotypic changes from other background-related 
variations (Soldner et al 2011; Reinhardt et al 2013). 
However, to date, no isogenic cardiac disease iPSC lines 
have been reported, and neither has the disease phenotype 
been investigated in parallel in both hiPSCs and hESCs 
harbouring the same mutation. 

In this study, we sought to generate and characterize a 
panel of control and LQT2-related human pluripotent stem 
cell (hPSC) lines. In order to investigate the specific effects of 
the C.A2987T (N996I) KCNH2 mutation without the con- 
founding elements that could result from individual genetic 
background variability, we studied this mutation under two 
genetically defined conditions. First, by genetically correcting 
the KCNH2 mutation in the LQT2-hiPSCs derived from a 
patient, and second, by introducing the same C.A2987T 
(N996I) mutation into an NKXZ.S^^^^^"^ hESC reporter 
line that allows the cardiac cells to be selected from a 
mixed differentiated population (ElHott et al 2011). 
Cardiomyocytes (CMs) differentiated from these two pairs 
of isogenic hPSCs demonstrated that the N996I KCNH2 
mutation causes an ~ 30-40% Jki reduction with 
consequential action potential (AP) duration (APD) 
prolongation, and indicated a trafficking defect as the 
associated pathological mechanism responsible for the 
disease electrophysiological phenotype. 

Results 

Generation of hiPSCs witli C.A2987T (N996I) KCNH2 
mutation 

Dermal fibroblasts were obtained from a 38-year-old woman 
with a diagnosis of famihal type-2 LQTS (Supplementary 
Figure SI A and B). Genetic screening indicated a heterozy- 
gous C.A2987T mutation in exon 13 of the KCNH2 gene 
(Figure lA), resulting in an asparagine to isoleucine substitu- 
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tion at position 996 (N996I) within the intracellular 
C-terminal region of the encoded HERG potassium channel 
(Figure IB). We generated hiPSCs (LQT2-hiPSCs^'^'^^^) from 
the primary skin fibroblasts by retroviral transduction of four 
reprogramming genes (0Cr4, S0X2, KLF4, and MYC) 
(Takahashi et al 2007a) (Figure IC), and three of the result- 
ing clones underwent further characterization. These clones 
expressed the hESC markers SSEA4 and NANOG, as detected 
by immunofluorescence staining (Figure ID). Additionally, 
qRT-PCR indicated that the hiPSC lines had reactivated the 
endogeneous pluripotency genes 0CT4, 80X2, LEFTYA, 
NANOG, REXl , and TDGFl , and silenced the retroviral trans- 
genes (Supplementary Figure S2 A and B) . Finally, differentia- 
tion of the hiPSCs as embryoid bodies demonstrated the 
capacity of all hues to generate derivatives of the three 
embryonic germ layers in vitro (Supplementary Figure S3). 
On the basis of these findings, one hiPSC clone was selected 
for targeted genetic correction. 

Targeted gene correction in LQT2-liiPSCs!^^^^' and 
targeted gene mutation in NKXZ.S^^'"''^'^ liESCs 

The study design is illustrated in Figure 2A. A conventional 
homologous recombination strategy was performed for either 
targeted gene correction in LQT2-hiPSCs^^^^^ or targeted gene 
mutation in the NKX2.5^^^^^'^ hESCs. We constructed two 
targeting vectors that included either a wild-type A or a 
mutated T nucleotide at KCNH2 c.2987 {KCNH2-A-loxP-pGK- 
Neo-loxP and KCNH2-T-loxP-pGK-Neo-loxP, respectively) 
(Figure 2B). Targeting of either the mutant allele (hiPSCs) 
or one of the two wild-type alleles (hESCs) is expected to 
result in a single base pair change that either repairs (LQT2- 
hiPSCs^°T or introduces (NKX2.5^^^^/"^ hESCs^'^'^^^) the 
N996I HERG mutation. The coding regions and the exon- 
intron boundaries of KCNH2 in the hiPSCs and hESCs were 
sequenced to exclude the presence of any other variants. 

The linearized targeting constructs were electroporated 
into the hiPSCs and hESCs and G418-resistant colonies iso- 
lated. Targeted clones were identified using a PCR-based 
screen that amplified the novel junction fragment between 
the genomic DNA and the integrated targeting vector gener- 
ated following homologous recombination (Figure 2B and C) . 
These PGR products were sequenced to confirm homologous 
recombination and to determine whether the mutation was 
corrected or introduced. We obtained one corrected 
LQT2-hiPSC clone (LQT2-hiPSCs'^°T and one mutated 
iVKX2.5"^^^/^ hESC clone (LQT2-hESCs^''''^^) (Figure 2D). 
The G418-resistance cassette was excised and the cell line 
cloned. Similar to the parental cell lines, the LQT2-hiPSCs''°^^ 
and NKX2.5'^^^^'^ hESCs^'^'^^^ expressed markers of undiffer- 
entiated hPSCs (Supplementary Figure S4A and B), could 
differentiate into the cellular derivatives of the three embryo- 
nic germ layers in vitro (Supplementary Figure S5), passed 
the PluriTest with a high 'pluripotency score' and a low 
'novelty score', indicating that they resemble normal hPSCs 
(Supplementary Figure S6A), and retained a normal karyotype 
(Supplementary Figure S6B and C). 

TAie C.A2987T (N996I) KCNH2 mutation causes a long-QT 
ptienotype in tiPSC-derived CMs 

The LQT2-hiPSCs"°'' and NKX2.5'^^^^'^ hESCs^''''^^ together 
with their parental cell lines, were differentiated to CMs, with 
spontaneously contracting foci first observed at 10-12 days of 

©2013 European Molecular Biology Organization 



Isogenic pairs of LQT2 pluripotent stem cells 

M Bellin et al 



KCNH2 



Wild-type reference 

CDS 2987A 

Protein sequence N996 

LQT2 patient 11-1 

CDS 2987A>T 

Protein sequence N996I 



E13 



B 






MA 


T C C A i 
S I 

AA/V\/ 


\ C A T T 
J I 

^aaAA 


T C C A ^ 
S [] 


J C A T T 



NH2 



HERG protein 




COOH 



LQT2-hiPSCsN996i 



LQT2-hiPSCsN996i 




Figure 1 Generation of hiPSCs from a patient with type-2 long-QT syndrome. (A) Genetic screening in the patient revealed the heterozygous 
single-nucleotide mutation A^T in exon 13 of the KCNH2 gene, in position 2987 of the coding sequence (CDS) (c.A2987T, NM_000238.3), 
resulting in the substitution of an asparagine with an isoleucine at position 996 of the protein (N996I, NP_000229.1). (B) The N996I mutation 
(red dot) is located in the C-terminal of the HERG protein, which is made of six trans-membrane domains (S1-S6), an amino (NH2) domain, a 
carboxyl (COOH) domain, and a pore (P) region. (C) Example of a hiPSC colony harbouring the C.A2987T (N9961) KCNH2 mutation (LQT2- 
hiPSCs^^^^^). Scale bar: 400 |im. (D) Immunofluorescence analysis of pluripotency markers SSEA4 (green) and NANOG (red) in a 
representative LQT2-hiPSC^^^^^ clone, with nuclear staining (DNA, blue). The image on the right is a magnification of the area framed in 
the left image. Scale bars: 100 iim (left image); 50 ^im (right image). 



differentiation. These were mechanically micro-dissected and 
allow^ed to mature further (Moretti et aU 2010). Dissociated 
cells displayed characteristic sarcomeric structures that w^ere 
positive for troponin I (TNNI) and alpha-actinin (Figure 3 A) , 
and continued to spontaneously contract. Quantitative RT- 
PCR in cardiac cells, corresponding to either micro-dissected 
beating areas (hiPSCs) or NKX2.5 eGFP+ cell populations 
(hESCs), revealed upregulation of both the KCNH2-la and -lb 
transcript variants, as w^ell as other ion channels essential for 
AP generation, compared to their undifferentiated counter- 
parts (Figure 3B). 

To confirm that the loxP sequence that remained in the 
KCNH2 locus following excision of the selection cassette did 
not interfere with HERG expression, we took advantage of 
another LQT2-hiPSC clone that was targeted and contained 
the residual loxP sequence, but was not corrected (LQT2- 
hiPSCs^''''^^/^°^P-"°^''°^). Analysing the response of the 
CMs to the selective Jkf channel blocker E-4031 using 
multielectrode arrays revealed that, consistent with earlier 
reports (Matsa et al, 2011), increasing concentrations of 
E-4031 caused prolongation of field potential duration 
(FPD) in all analysed LQT2-hiPSC lines (original LQT2- 
hiPSCs^'^'^'^ LQT2-hiPSCs^'^'^'^/^°^P-^°^^^°\ and LQT2- 
hiPSCs''°^ (Supplementary Figure S7A). Importantly, the 
field potential response did not differ between LQT2- 
j^.p5(.N996i LQT2-hiPSC^'^'^^/^°^P-^°^^^°^ CMS and was pro- 
longed when compared to LQT2-hiPSC^°'"'" CMs, indicating 
not only that the N996I mutation conferred increased sensi- 
tivity to Jki blockade but also that the remaining loxP 
sequence did not appear to interfere with their electrophy- 
siological phenotype. We therefore used the LQT2- 



hiPSCs^''''^^/^°^P""°^''°^ line for further electrophysiological 
characterization. Moreover, similar analysis in CMs from 
iVKX2.5^^^^/"^ hESCs and NKX2.5^^^^/^ hESCs^'^''^ showed 
that introduction of the N996I mutation in a genetically 
distinct background resulted in a significant, E-4031 -depen- 
dent FPD prolongation that was comparable to that measured 
in mutated LQT2-hiPSC CMs, suggesting an exclusive effect 
of this mutation on the observed phenotype (Supplementary 
Figure S7A-D). 

To further investigate whether the correction of the 
C.A2987T (N996I) mutation resulted in an increase in 
the repolarizing potassium current specifically transduced 
by the HERG channel, we recorded Jkf in individual hiPSC- 
derived CMs. Typical examples of current traces recorded 
before and after addition of 1 jiM of E-4031, in mutated and 
corrected LQT2-hiPSC-derived CMs are shown in Figure 4A 
and B, where % is defined as E-4031 -sensitive current. 
Average % density in LQT2-hiPSC^^^^^ CMs was significantly 
smaller than in the LQT2-hiPSC'^°^^ CMs (Figure 4C and D). At 
a membrane potential of — 20 mV, both Jki density (measured 
at the end of the test pulse) and peak-tail % density were 
decreased by 33 % . To confirm whether the C.A2987T (N996I) 
mutation is the cause of the LQT2 phenotype observed in 
hPSC-derived myocytes, we measured Jkf in CMs derived 
from NKX2.S'^^^'" hESCs and NKX2.5"^^^/^ hESCs^''''^^ that 
represent a different genetic background. Representative ex- 
amples of current traces recorded in wild-type and mutated 
hESC-CMs are shown in Figure 5A and B. Average Iki density 
was significantly reduced in the mutated hESC-CMs com- 
pared with the wild-type cells (Figure 5C and D). At a 
membrane potential of -20mV, Jki density (measured at 
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Figure 2 Gene targeting by homologous recombination in LQT2-hiPSCs^'^'^^^ and in NKX2.5^^^^^^ hESCs. (A) Schematic showing the project 
rationale. Patient-specific LQT2-hiPSCs harbouring the N9961 mutation were corrected and NKX2.5^^^^^^ hESCs were mutated by gene 
targeting. Parental and genetically modified hPSC lines were differentiated into CMs and their electrophysiological phenotypes were analysed 
and compared. (B) The strategy for precise genomic modification of KCNH2. Top line, structure of the KCNH2 locus. Numbered black boxes 
indicate exons 6-15. Exon 13 is mutated (red T) in LQT2-hiPSCs^^^^^ and wild type (black A) in NKX2.5''^^^^'" hESCs. The gene targeting vector 
for correcting the mutation in LQT2-hiPSCs has the wild-type adenine nucleotide, whereas the gene targeting vector for introducing the 
mutation in NKX2.5^'^^^^^ hESCs has the mutated thymine nucleotide. NeoR, the PGK-Neo cassette encoding G418 resistance flanked by loxP 
sequences (grey triangles), was inserted in the reverse direction. PGR primers (a, b) and (c, d) were used to identify the targeted clone. (C) PGR 
analysis using these primers generated specific bands of 4.8 kb (5' homology arm) and 7.3 kb (3' homology arm) from targeted clones (LQT2- 
hiPSCs^^°^ and iVXX2.5^^^^/"' hESGs^^°^). (D) Sequence analysis of PCR-amplified genomic DNA showing correction of the C.A2987T mutation 
in the LQT2-hiPSGs''°'' line and mutation in the NKX2.S''^^^^'" hESCs^'^'^^^ line. The wild-type reference sequence (Ref) is shown in the top line. 
Source data for this figure is available on the online supplementary information page. 



the end of the test pulse) and peak-tail Jki density were 
decreased by 43 and 45%, respectively, corroborating the 
results obtained in the patient-specific hiPSC model. 
Importantly, while Iki densities were affected by the muta- 
tion, voltage dependence of activation parameters (half-max- 
imal voltage [Vi/i) and slope factor (/c)) were similar, both in 
mutated and corrected LQT2-hiPSC-derived CMs and in wild- 
type and mutated hESC-derived CMs (Figure 6A and B). 
Furthermore, the time constant of Jkf activation (x) 
(Figure 6C and D), as well as the fast (if) and the slow (xg) 
deactivation time constants were all unchanged when com- 
paring wild-type and mutated cells (Figure 6E and F) . Taken 
together, these findings suggest that the N996I HERO muta- 
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tion is the only relevant mutation that influences Jkf current 
in the heterozygous LQT2 CMs. 

To evaluate the impact of the Jki change on the APD, APs 
were recorded in single spontaneously contracting cells 
(Figure 7) . Representative APs (paced at 1 Hz) from mutated 
and corrected hiPSC-derived CMs, and from wild-type and 
mutated hESC-derived CMs are shown in Figure 7A and C, 
respectively. The APD at 50 and 90% of repolarization 
(APD50 and APD90, respectively) was significantly reduced 
in the LQT2-hiPSC''°'' CMs compared with their mutated 
counterparts (Figure 7B). In the hESC model, only the 
APD50 of NKX2.5'^^^'^ hESC^^^^^ CMs was significantly 
prolonged compared with the NKX2.5^^^^^'^ hESC CMs 
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Figure 3 Differentiation of hPSCs into the cardiac lineage. (A) Confocal immunofluorescence images of cardiac sarcomeric proteins TNNI 
(green) and a-actinin (red) in human CMs generated from LQT2-hiPSCs^^^^^ LQT2-hiPSCs^°'', NKX2.5^^^^/^ hESCs, and Ni^X2.5^^^^/^ 
hESCs^^^^^ Nuclei are stained in blue. Bottom panels are a magnification of the area framed in the upper images. Top panels, scale bar: 25 |im; 
bottom panels, scale bar: 10 |im. (B) Transcriptional profile of human CMs generated from mutated and corrected hiPSCs (LQT2-hiPSCs^^^^^ 
and LQT2-hiPSCs^°'', respectively) and from wild-type and mutated hESCs (Ni^X2.5^^^^/^ hESCs and NKXZ.S''^^^^'" hESCs^^'^^^ respectively). 
Undifferentiated cells from each hPSC line are also shown (Undiff.) . Quantitative RT-PCR analysis was performed on the cardiac troponin gene 
(TNNT2), to show enrichment for the cardiomyocyte population, on the HERG channel gene [KCNH2], and on other key genes encoding for ion 
channels involved in the generation of the action potential in cardiac cells. All values are normalized to GAPDH and are relative to 
undifferentiated NKX2.5^^^^^ hESCs. Raw minimum (min) and raw maximum (max) values were taken as a reference for heatmap 
representation. 
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Figure 4 densities in mutated and corrected LQT2-hiPSC-derived CMs. (A) Representative current traces elicited upon 4 s depolarizing 
voltage steps to - 30, - 20, and - 10 mV from a holding potential of - 40 mV, before and after the application of 1 |iM E-4031. Inset: voltage 
protocol. (B) Typical tail currents measured after the depolarizing step to - 10 mV in corrected (black) and mutated (red) LQT2-hiPSC-CMs, 
showing a bi-exponential decay. (C) Average current-voltage [I-V] relationships for /ki^ measured at the end of the test pulses, in mutated (red) 
and corrected (black) LQT2-hiPSC-derived CMs. Inset: voltage protocol. * indicates statistical significance (P= 0.046, two-way rmANOVA; 
Holm-Sidak test post hoc analysis: -20mV: P=0.015, - lOmV: P = 0.007, OmV: P = 0.018, lOmV: P = 0.011, 20mV: P = 0.013). (D) Average 
7- y relationships for peak tail currents in mutated (red) and corrected (black) LQT2-hiPSC-derived CMs. * indicates statistical significance 
(P = 0.039, two-way rmANOVA; Holm-Sidak iesi post-hoc analysis: - 20 mV: P = 0.008, - 10 mV: P = 0.004, 0 mV: P = 0.005, 10 mV: P = 0.004, 
20 mV: P = 0.006). Inset: voltage protocol. 
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(Figure 7D). In contrast, in both the hiPSC and hESC para- 
digms, no differences in upstroke velocity [V^^^], AP amph- 
tude (APA), and maximal diastolic potential (MDP) were 
detected between the mutated and corrected patient hues, 
and between the wild-type and mutated hESC groups 
(Figure 7B and D, respectively). 

The C.A2987T (N996I) KCNH2 mutation causes a 
trafficking defect in ttie HERG ciiannel 

To investigate the functional consequences of the N996I- 
HERG mutation further, we examined the cellular distribution 
of the HERG channel in hiPSC- and hESC-derived CMs. 
Immunocytochemistry in both hPSC-CMs revealed a distribu- 
tion of the HERG channel over the cell surface (Figure 8 A) , as 
well as in the endoplasmic reticulum (ER) and in the Golgi 
apparatus (Supplementary Figures S8 and S9). However, it 
was not possible to determine accurately by immunofluores- 
cence whether CMs generated from the LQT2-hPSCs^^^^^ 
models presented any differences in subcellular localization 
of HERG channels when compared to their wild-type or 
corrected counterparts. Therefore, we exploited the fact that 
we had introduced the C.A2987T (N996I) KCNH2 mutation 
into an NKXZ.S^^^^^"^ hESC line in which GFP reports the 
expression of the cardiac associated transcription factor 
NKX2.5 (Elliott et al, 2011). This allowed us to purify a 
population enriched for CMs by flow cytometry (Figure 8B), 
and quantify HERG channel levels within this fraction by 
western blot analysis. Two HERG protein bands were identi- 
fied consistent with core glycosylated (135 kDa, immature) 
and more complex glycosylated (155 kDa, mature) isoforms 



(Figure 8C), with the latter present predominantly in the 
NKX2.5-eGFP+ fraction. In the iVKX2.5-eGFP+ hESC^^^^^ 
cells, the 155-kDa protein band, representing the form trans- 
ported to the cell membrane through the Golgi, was reduced 
by ~ 2-fold compared to the wild-type cells, while the 
135-kDa band, which corresponds to the protein located in 
the ER, was unaffected (Figure 8D). This resulted in a 
significant decrease in HERG trafficking efficiency, calculated 
as a ratio of the fully glycosylated 155 kDa band over total 
HERG protein (155 kDa band + 135 kDa band), of -40% 
(Supplementary Figure SlOA). Similar analysis in LQT2- 
hiPSC"°'' and LQT2-hiPSC^'^'^^^ CMs confirmed a specific 
reduction in the complex glycosylated HERG protein in the 
mutated cells compared to the corrected counterparts 
(Figure 8E and F; Supplementary Figure SlOB). Thus, the 
N996I mutation appears to interfere with the maturation and 
trafficking of the channel to the membrane. 

To further investigate the pathological mechanism of the 
N996I mutation, we tested whether this mutation resulted in 
a misfolded HERG protein by analysing the unfolded protein 
response (UPR) in CMs purified from NKXZ.S^^^^^"^ hESCs 
and NKX2.5^^^^/^ hESCs^'^'^^^ UPR is an ER stress pathway 
that increases the synthesis of chaperones, such as Calnexin 
and Calreticuhn, which in turn target misfolded proteins for 
proteasome degradation (Ellgaard and Helenius, 2003; Hetz, 
2012). Western blot for the Activating Transcription Factor 6 
(ATF6), a key regulator of transcriptional control in the 
mammalian UPR (Haze et al, 1999; Chen et al, 2002), 
revealed the presence of several forms of ATF6, including 
both the unprocessed protein embedded in the ER (90 kDa 
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Figure 6 /kp activation and deactivation properties in mutated and corrected LQT2-hiPSC-derived CMs and in wild-type and mutated hESC- 
derived CMs. (A, B) Average peak tail current normalized to the maximal current following repolarization to -40mV in mutated (red) and 
corrected (black) LQT2-hiPSC-derived CMs (A) and in wild-type (black) and mutated (red) hESC-derived CMs (B). Inset: voltage protocol. 
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band) and some smaller forms resulting from the cleavage of 
the full-length ATF6 during ER stress, with no detectable 
differences between wild-type and mutated MCX2.5-eGFP + 
cells (Figure 9A). Moreover, there was no evidence of an 
increase in the expression of the chaperones Calnexin 
and Calreticuhn in NKXZ.S^^^^^"^ hESC^'^'^^^ CMs (Figure 9B 
and C), indicating that the N996I HERG mutation does not 
appear to stimulate the UPR pathway. Nevertheless, to assess 
the involvement of proteasome-mediated and/or lysosome- 
mediated degradation of HERG channel further, we deter- 
mined the effects of the proteasome inhibitor lactacystin and 
the lysosome inhibitor leupeptin on the steady-state HERG 
protein levels in wild-type and mutated hESC-CMs. Following 
purification of the M<X2.5-eGFP^ cell population by flow 
cytometry, cells were treated for 24 h with either lactacystin 
(20|^M) or leupeptin (100 |iM) and cell lysates were then 
immunoblotted with anti-HERG antibody (Figure 9D). 
Treatment with the proteasome inhibitor lactacystin augmen- 
ted the level of the mature, highly glycosylated form of HERG 
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in NKXZ.S'^^^'^ hESC^''''^^ eGFP+ cells, resulting in an ~ 1.4 
fold increase in channel trafficking efficiency. No increase 
was observed in the wild-type counterparts or after treatment 
with the lysosome inhibitor leupeptin (Figure 9D and E), 
suggesting a role of proteasomes in the degradation of N996I 
HERG protein. 

Discussion 

The ability to reprogramme somatic cells from any individual 
of choice has created unprecedented opportunities to study 
disease development in vitro and create novel platforms for 
drug screening and discovery. However, one crucial limita- 
tion for disease modelling in hiPSCs has been that, unlike 
experimental animal models, humans are genetically hetero- 
geneous (Goldstein, 2009; McKernan et al, 2009), and 
genetically matched controls are necessary to distinguish 
disease-relevant phenotypic changes from normal 
background-related variations due to polymorphisms in 
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other genes. Here, we generated a panel of isogenic control 
and cardiac disease lines and investigated the disease 
phenotype in parallel in both hiPSCs and hESCs harbouring 
the same N996I KCNH2 mutation. Although this mutation has 
been described as a novel LQT2 variant and included in 
studies exploring possible genotype-phenotype correlations 
(Khositseth et al 2004; Tester et al 2005; Tan et al 2006), its 
mechanism of action has remained unclear. Here, we 
demonstrated that this specific mutation is the cause of the 
LQT2 phenotype observed in the hPSC-derived CMs. 
Furthermore, we investigated the biophysical disease 
mechanism in functionally relevant human cardiac cells. 
Our findings demonstrated that the N996I KCNH2 mutation 
caused a similar ( ~ 30-40 % ) Ikf reduction in both mutated 
hiPSC- and hESC-derived CMs. In particular, because a 
decrease of ^50% in the cardiac repolarizing potassium 
current indicates haploinsufficiency (Moss et al 2007; 
Vandenberg et al 2012), our results suggested this as the 
mechanism underlying Jkf reduction, rather than a dominant- 
negative effect of mutant subunits on the wild-type subunits. 
Furthermore, the reduction in the fully glycosylated mature 
155 kDa HERO protein in the mutated hPSC-CMs and its 
partial rescue after treatment with the proteasome inhibitor 
lactacystin suggest that the mutation might induce, to some 
extent, protein clearance through the proteasomes, thus 
resulting in defective channel trafficking (Vandenberg et al 
2012). Further analysis of HERO protein stability and 
internahzation is required to unravel the additional 
mechanisms through which mutated N996I HERO channels 
may be processed. It is noteworthy that, in contrast to 
previous reports examining trafficking-deficient HERO 
mutants in heterologous systems (Smith et al 2011; Quo 
et al 2012; Wang et al 2012), we did not observe any 
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accumulation of the core-glycosylated form of HERO in the 
ER nor an activation of the UPR in the mutated hPSC-derived 
CMs. This might reflect either the specific effect of this 
particular mutation or the different cell systems used to 
investigate the pathophysiological mechanism. In 
overexpression systems, the levels of both wild-type and 
mutated HERO proteins are far from physiological; this 
could alter molecular stoichiometry and cause a more 
pronounced stimulation of ER stress and clearance 
pathways than in the native conditions of a CM. 

It is widely recognized that hPSC-derived CMs are imma- 
ture, with electrical properties, gene expression and contrac- 
tion forces only equivalent to those of fetal or neonatal 
human myocytes (de Boer et al 2010; Hoekstra et al 2012). 
Nevertheless, our results show that these cells are capable of 
capturing specific traits of an electrical disease of the heart, 
such as LQT2, corroborating earlier findings (Itzhaki et al 
2011; Matsa et al 2011). Of note, AP parameters (especially 
upstroke velocity) suggested that hESC-derived CMs were 
electrophysiologically less mature than those derived from 
the hiPSCs used here. While this did not preclude detection of 
Jkf reduction as a specific effect of the N996I-HERG mutation, 
the electrophysiological immaturity of the hESC-CMs might 
have influenced the impact of Jxr reduction on APD, where 
we observed a significant prolongation of APD50 but not of 
both APD50 and APD90 as observed in hiPSC-CMs. In this 
regard, the longer plateau phase of the AP in hiPSC-CMs 
compared to hESC-CMs would allow a higher number of 
HERO channels to inactivate, resulting in a different 
contribution of Jki to AP repolarization (Rudy, 2008) . These 
differences could also be due to inherent variability among 
individual hiPSC and hESC lines (Narsinh et al 2011; 
Mummery et al 2012) but do not detract from the 
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conclusion since isogenic pairs were always compared. In 
addition, while the Jkf traces recorded in our lines are in 
agreement with % previously reported in hPSC-derived CMs 
(Moretti et al 2010; Itzhaki et al 2011; Ma et al 2011; Lahti 
et al 2012), the current density measured in our hESC-CMs 
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was higher (3.5 pA/pF) than those described in the literature 
for both hiPSC-derived (0.55-1.9 p/pF) and native human 
(0.25-0.6 pA/pF) CMs (Hoekstra et al 2012). The molecular 
basis of this discrepancy requires further independent 
investigation, but a larger Jkf could impact on the AP 
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characteristics, for example, contributing to a faster 
repolarization resulting in a shorter AP. Indeed, while the 
APD values measured in the LQT2-hiPSC''°^^ CMs fall into the 
range reported for hPSC-derived CMs (Hoekstra et al, 2012), 
in the NKX2.5^^^^/^ hESC CMs the APD is rather small. 

Another aspect of our work is noteworthy: we measured 
very modest APD prolongation in both LQT2-hiPSC^'^'^^^ and 
NKXZ.S^^^^^"^ hESC^'^'^^^ CMs when compared to their wild- 
type counterparts, with values shorter than those reported in 
other LQT2 hiPSC models (Itzhaki et al, 2011; Matsa et al, 
2011; Lahti et al, 2012), and did not observe any early-after 
depolarization in the mutated cells. These differences are 
hkely due to the specific N996I mutation, which results in a 
less severe LQT2 phenotype than those caused by other 
trafficking-deficient mutations (Moss et al, 2002). However, 
it should be noted that in all previous studies using LQT2 
hiPSC models the control group consisted of wild-type hiPSCs 
derived from a healthy individual unrelated to the LQT2 
patient, therefore over-estimation of differences cannot be 
excluded under these circumstances. The advantage of using 
isogenic hues is indeed highhghted by our results, where 
even subtle differences could be appreciated. 

Most importantly, our data provide evidence that the 
pathogenesis of the N996I-HERG mutation can be modelled 
exclusively in CMs generated from hPSCs, without the need 
for heterologous and/or overexpression systems. In particu- 
lar, we detected and analysed physiological levels of HERG 
protein (by immunostaining and western blot), as well 
as molecular chaperones and transcriptional regulators of 



the UPR, in functionally relevant human cardiac cells. 
Moreover, the use of an NKXZ.S^^^^^"^ hESC reporter line 
allowed selection of the cardiac cell population from other 
cell types, such as neurons, endocrine cells, and undiffer- 
entiated proliferating cells, that also express KCNH2 (Rosati 
et al, 2000; Pillozzi et al, 2002; Huffaker et al, 2009; 
Vandenberg et al, 2012). This highUghts the particular 
value of combining patient-specific iPSC technology and 
genetic modification tools. 

Since our study compared genetically matched lines differ- 
ing uniquely in one LQT2 point mutation, a clear genotype- 
phenotype relationship could be ascertained. This is of 
particular value in avoiding over- or under-estimation of 
the impact of disease-related changes as a result of inadver- 
tent bias when choosing a random control (even if age and 
gender matched with the same reprogramming and differen- 
tiation method). Therefore, this approach is a powerful 
means of identifying links between specific genotypes and 
cardiac disease predisposition, as also recently demonstrated 
for neurodegenerative disease (Reinhardt et al, 2013). 

Finally, targeted gene correction of KCNH2 in hiPSCs rescued 
Jki density in hiPSC-CMs and normaUzed AP values. Although 
precise gene modification of human CMs in situ is not an 
immediate option for regenerative medicine or for treating 
LQTS, gene targeting as carried out here provides essential 
tools for deciphering molecular mechanisms, studying genetic- 
phenotypic interaction and perhaps, in the future, analysing 
the potential disease-causing SNPs identified by genome-wide 
association studies in humans. Other rescue approaches, such 
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as RNA interference, may only result in partial knock down of 
the mutated mRNA (Matsa et al, 2013], while genetic 
correction should restore the normal genetic phenotype in 
monogenic diseases. This can then demonstrate whether the 
specific genetic lesion is pathogenic or whether there are likely 
to be other modifier loci involved. Considering the increasing 
awareness of the broad influence of genetic background 
(Montgomery et al, 2010] and the biological differences 
between hPSC lines (Boulting et al, 2011], our experimental 
system may overcome some of the shortcomings of 
conventional hiPSC approaches in identifying disease-related 
phenotypes. As more LQTS-hiPSC lines are investigated using 
this approach, the molecular and genetic reasons for variability 
of expression and for incomplete penetrance seen in this, and 
other channelopathies, may be elucidated. 

In summary, we have demonstrated that genetic correction 
of the N996I KCNH2 mutation associated with LQT2 restores 
Jki density and normalizes APD in patient-specific LQT2- 
hiPSC-CMs. The same mutation introduced into a different 
genetic background (hESC] results in a very similar electro- 
physiological phenotype. Furthermore, the reduction in the 
mature complex-glycosylated 155 kDa protein band in the 
mutated hPSC-CMs suggests that, due to the mutation, a 
smaller number of channels are present at the cell membrane. 
This could be partially due to proteasome targeting for 
degradation, but might not be the only mechanism for this 
N996I-specific mutation. Further studies are required to 
assess more in detail the exact contribution of proteasome 
clearance pathway, protein channel stabiUty, and lysosome 
degradation. Our study demonstrates that isogenic pairs of 
hPSCs can be used (1] to prove the authenticity of the 
genotype-phenotype correlation and (2] to unravel the 
pathophysiological mechanism seen in a genetically inherited 
cardiac disease. As such, this approach represents a robust 
strategy to study not only cardiac disease mechanisms but 
also other genetic disorders. 

Materials and methods 

Clinical history and genetic phenotype 

A 38-year-old Caucasian woman presented with a prolonged QT 
interval (QT interval corrected for heart rate (QTc] 617 ms] at the 
ECG and was diagnosed with type-2 LQTS. Genetic screening 
indicated a heterozygous mutation C.A2987T in exon 13 of the 
KCNH2 gene, resulting in the N996I missense mutation. The patient 
has thus far been asymptomatic. She is now treated with p-blockers. 



Generation of patient-specific hiPSCs, differentiation into 
CMs, and drug treatment 

We recruited a 38-year-old Caucasian female LQT2 patient for 
dermal biopsy after obtaining written informed consent. 
Reprogramming of primary skin fibroblasts was performed as 
described previously (Moretti et al, 2010; Jung et al, 2012]. 
Briefly, fibroblasts were infected with retroviruses encoding 0CT4, 
S0X2, KLF4, and MYC and cultured on mouse embryonic feeder 
(MEF] cells until hiPSC colonies could be picked. Both hiPSCs and 
hESCs were maintained in culture using standard procedures 
(Takahashi et al, 2007a] and were differentiated using standard 
cardiogenic protocols (Mummery et al, 2003; Elliott et al, 2011]. 
Spontaneously contracting areas were manually dissected and 
cultured further until days 20-60 of differentiation. Cells for physio- 
logical experiments were enzymatically dissociated into single cells 
(Moretti et al, 2010], plated on fibronectin-coated glass coverslips, 
and analysed within 7-15 days. 

For drug treatment, NKX2.5-eGFP^ cell populations were purified 
by flow cytometry and maintained in low attachment plates while 



treated for 24 h with 20 |iM lactacystin (Sigma] or 100 |iM leupeptin 
(Sigma]; cells were then lysated and subjected to immunoblotting. 

Genomic sequencing 

Genomic DNA was isolated from cultured cells using the Centra 
PureCene Cell Kit (Qiagen] . To confirm the presence or absence of 
the C.A2987T KCNH2 mutation in the hPSCs, and to exclude the 
presence of other variants in the KCNH2 locus, the coding region 
and the exon-intron boundaries were PCR amplified from genomic 
DNA (Supplementary Table SI]. The PCR products were purified 
using the QIAquick PCR Purification kit (Qiagen] and sequenced. 

Generation and identification of targeted hiPSCs and hESCs 

A BAC (RP11-10L20] containing the -33-kb human KCNH2 locus 
was modified stepwise by recombineering to generate the final 
targeting vector (Fu et al, 2010]. First, a /oxP-flanked positive 
selection cassette [loxP-PGK-neo-loxP] comprising a mammalian 
promoter [pGK], a bacterial promoter IgbS), and a G418/ 
kanamycin-resistance gene [neo] was amplified from an R6K 
plasmid and inserted 0.9 kb downstream of the mutation site. A 
15.6-kb fragment, including a 6.6-kb 5'-homology arm, the inserted 
G418-resistance cassette, and a 4.3-kb 3'-homology arm, was then 
sub-cloned into a minimal vector (Gene Bridges], generating the 
wild-type KCNH2 targeting vector [KCNH2-A-loxP-pGK-Neo-loxP] 
used to correct the mutation in the LQT2-hiPSCs^^^^\ For hESC 
targeting, the C.A2987T mutation was introduced into this targeting 
vector using the QuickChange XL Site-Directed Mutagenesis Kit 
(Stratagene] [KCNH2-T-loxP-pGK-Neo-loxP) . The two vectors, wild 
type and mutated, were linearized with the restriction enzyme Notl 
before electroporation (Costa et al, 2007]. Targeted clones were 
identified using a PCR-based screening strategy. Homologous 
integration of both the 5' and 3' homology arms, as well as 
correction or introduction of the C.A2987T KCNH2 mutation in the 
resulting clones were confirmed by sequencing. The /oxP-flanked 
G418-resistance cassette was excised using Cre recombinase from 
both a corrected- and non-corrected hiPSC clone, and a mutated 
hESC clone, and the resulting lines subcloned as described 
previously (Davis et al, 2008]. Clonal lines were screened for the 
loss of the neomycin-resistance cassette and for the absence of 
the Cre recombinase expression plasmid by PCR (Supplementary 
Table SI]. These PCR products were also sequenced to confirm 
excision of the neomycin-resistance cassette. Karyotype analysis 
was performed using COBRA-FISH as described elsewhere (Szuhai 
and Tanke, 2006]. Twenty metaphase spreads for each sample were 
analysed. 

Immunofluorescence and western blot analysis 

Cells were fixed in 4% paraformaldehyde, permeabilised with 
phosphate buffer saline (PBS]/0.1% Triton X-100 (Sigma-Aldrich] , 
and blocked with 10% FCS (Life Technologies]. Samples were 
incubated overnight at 4°C with primary antibodies specific for 
the following: NANOG (goat polyclonal, R&D Systems], SSEA4 and 
TNNI (mouse monoclonal and rabbit polyclonal, respectively, both 
from Santa Cruz Biotechnology], a-ACTININ and HERG (mouse 
monoclonal and rabbit polyclonal, respectively, both from Sigma- 
Aldrich], PDI (mouse monoclonal. Abeam], and Golgin-97 (mouse 
monoclonal. Molecular Probes]. Primary antibodies were detected 
with either Cy3- or Alexa-Fluor 488-conjugated antibodies. Nuclei 
were visualized with DAPI (Invitrogen] and F-actin with Phalloidin- 
Alexa-Fluor-594-conjugate (Invitrogen]. Images were captured 
using either a Leica DMI6000-AF6000 fluorescence microscope or 
a Leica SP5 confocal laser-scanning microscope (both from Leica 
Microsystems]. 

Western blotting on whole-cell lysate of either eCFP ^ and eCFP ~ 
cells purified from differentiated iVXX2.5^^^^/"' hESCs and 
NKX2.5''^^^^'^ hESCs^'^'^^^ or of manually microdissected beating 
areas from differentiated LQT2-hiPSCs^^^^^ and LQT2-hiPSCs^°''" 
was performed with standard protocols (Kurien and Scofield, 
2006] using 40-80 |ig proteins and the following primary anti- 
bodies: HERG (mouse monoclonal, Enzo Life Sciences], ACTIN 
(mouse monoclonal, Millipore], ATF6 (mouse monoclonal, 
ActiveMotif] , Calnexin and Calreticulin (both rabbit polyclonal. 
Abeam]. Western blots were quantified using the ImageJ software 
(http://rsb.info.nih.gov/ij/]. Whole film images are provided in 
Supplementary data. 
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Gene expression analysis 

For qRT-PCR, total RNA was purified using the RNeasy Mini Kit 
(Qiagen) according to the manufacturer's protocol. RNA samples 
were reverse transcribed using the Superscript 11 First-Strand 
Synthesis kit (Invitrogen) . Gene expression was assessed by qRT- 
PCR as described previously (Moretti et al, 2010). Gene expression 
levels were normalized to GAPDH and represented in a heatmap 
using the GENE-E software (http://www.broadinstitute.org/cancer/ 
software/GENE-E/index.html). Primer sequences are provided in 
Supplementary Table SI. 

For whole-genome microarray, total RNA was isolated from 
undifferentiated hPSCs using the NucleoSpin miRNA kit 
(Macherey-Nagel) and hybridized on lllumina HT12v4 microarrays. 
The raw microarray data were analysed with the PluriTest algorithm 
(http://www.pluritest.org; Muller et al, 2011). 

Elect rophi ysiological characterization 

Data acquisition and analysis. Recordings were performed on 
single CMs 7-15 days after cell dissociation. For Iki measurements, 
CMs were identified on the basis of their typical morphology, while 
for AP measurements, spontaneously contracting cells that could be 
paced at the stimulation frequency of 1 Hz were selected. Data were 
collected from at least three independent differentiations per line. 
/ki and APs were recorded with the ruptured and the perforated 
patch-clamp technique, respectively, using an Axopatch 200B am- 
plifier (Molecular Devices). Voltage control and data acquisition of 
/ki and APs were performed with pClamplO.l (Axon Instruments). 
Membrane currents were analysed with Clampfit 10.1 (Axon 
Instruments), while analysis of APs was performed with the cus- 
tom-made software. Patch pipettes had a tip resistance of ~2MQ 
and series resistance (Rs) was compensated by 80%. /kp and APs 
were filtered (2 and 5 kHz, respectively), and digitized (2 and 
40kHz, respectively). Potentials were corrected for liquid junction 
potentials (Barry and Lynch, 1991), which were calculated using 
Clampex 10.1 (Axon Instruments) and were equal to 15 and 17mV 
for AP and /kp measurements, respectively. 

Voltage-clamp experiments. iKr was recorded at 37°C using a 
pipette solution containing (mM): K-gluconate 125, KCl 20, K2- 
ATP 5, HEPES 10, EGTA 10; pH 7.2 (KOH). The bath solution was 
Tyrode's solution containing (mM): NaCl 140, KCl 5.4, CaCl2 1.8, 
MgCl2 1.0, glucose 5.5, HEPES 5; pH 7.4 (NaOH). Nifedipine (5 |iM) 
and JNJ-303 (1 |iM) were added to the external solution to block the 
L-type calcium current (/caJ and the slow component of the 
rectifier potassium current (7ks). respectively. Currents were elicited 
using 4 s hyper- and depolarizing pulses from a holding potential of 
- 40 mV, as indicated in the voltage protocols depicted in Figures 4 
and 5. The cycle length of the voltage clamp steps was 10 s. Iki was 
measured as a l-|iM E-4031 -sensitive current, by subtraction of the 
current recorded before and after E-4031 application. densities 
were calculated by dividing current amplitude (pA), measured 
either at the end of the test pulses (Figures 4C and 5C) or at the 
peak of the tail current (Figures 4D and 5D), by cell membrane 
capacitance (pF). Cell membrane capacitance was measured by 
dividing the decay time constant of the capacitive transient in 
response to 5mV depolarizing steps from -40mV, by the Rs. 
Average cell capacitance for both hiPSC- and hESC-CM groups are 
shown in Supplementary Figure Sll. 

Activation curves, determined from peak tail current normalized 
to the maximal value and plotted against the test pulse voltage, 
were fitted with Boltzmann equation (y = [1 + exp{(y- ^1/2)/ 
fc}]~^), where V1/2 is the half-maximal voltage of activation, and 
k is the slope factor. 

Time constants of activation (x) were obtained by fitting /ki 
during depolarizing pulses with a single-exponential function 
(y = Ao + Aiexp[ - t/x]). Finally, the slow (xg) and fast (xf) time 
constants of deactivation were determined by fitting /ki tail currents, 
at voltage steps positive to - 10 mV, with a bi-exponential function 
(y = Ao + Afexp[ - t/xf] + Asexp[ - t/Xg]) . 

Current-clamp experiments. APs were recorded at 37°C in Tyrode's 
solution. The pipette solution contained (mM): K-gluconate 125, 
KCl 20, NaCl 5, amphotericin-B 0.22, HEPES 10; pH 7.2 (KOH). APs 
were elicited at the stimulation frequency of IHz by 3 ms, 1.2 x 



threshold current pulses through the patch pipette. MDP, Vmax. APA, 
and APD50 and APD90 were analysed. Data from 10 consecutive APs 
were averaged. 

MEA electrophysiology. MEA standard measurements were per- 
formed in DMEM culture medium supplemented with 2% FCS as 
previously described (Braam et al, 2013). E-4031 was dissolved in 
DMSO at 10 mM and serial dilutions were made in culture medium. 

Statistical analysis 

The SigmaStat 3.5 software was used for statistical analysis. Results 
are expressed as mean±s.e.m. The software performed tests for 
normality and equality of variance on the data sets prior to 
application of a statistical test. Group comparisons were made 
using unpaired Student's Mest when both conditions were fulfilled 
or Mann- Whitney test when the conditions of normality and equal 
variance were not met. For repetitive measurements, two-way 
rmANOVA followed by a Holm-Sidak test for post-hoc analysis 
was used. P<0.05 was considered as statistically significant. 

Ethics statement 

The protocols for research involving human subjects and for stem- 
cell research were approved by the institutional review board and 
the committee charged with oversight of embryonic stem-cell 
research at the Technical University of Munich. Written informed 
consent was received from the participant prior to inclusion in the 
study. The studies on hESCs were performed exclusively in the 
Netherlands where the ethical committee for research at the LUMC 
provided approval for their use. 

Supplementary data 

Supplementary data are available at The EMBO Journal Online 
(http://www.embojournal.org) . 



Acknowledgements 

We thank C Hohnke (Technical University of Munich) for obtaining 
the skin biopsy; K Takashi and S Yamanaka (Kyoto University) 
for providing viral vectors through Addgene; F Stewart (Technical 
University of Dresden) for providing the R6K plasmid; MJM van der 
Burg, K Szuhai, and H Tanke (Leiden University Medical Center) for 
karyotyping analysis; AO Verkerk (University of Amsterdam), 
M Birket (Leiden University Medical Center), D Sinnecker 
(Technical University of Munich) and R Dirschinger (Technical 
University of Munich) for comments on the manuscript; S 
Bohringer (Leiden University Medical Center) for statistical 
assistance. We are especially thankful to the type-2 long-QT 
syndrome patient for tissue donation. This work was supported 
by grants from the European Research Council, ERG 261053 
(K-LL) and ERG 323182 (CLM); the German Research Foundation, 
Research Unit 923, Mo 2217/1-1 (AM), La 1238 3-1/4-1 (K-LL), 
and WE 1639/3-1 (AW); German Centre for Cardiovascular Research 
(K-LL and AM); ZoNMW Animal Alternatives, 114000101 (CLM, 
SC, and LT); Netherlands Proteomics Consortium, 050-040-250 
(CLM); Netherlands Institute of Regenerative Medicine (CLM and 
RD); EU Marie Curie FP7-people-2011-lEF programme, HPSCLQT 
29999 (MB) and a short-term EMBO fellowship, ASTF 387.00-2011 
(CD'A). 

Author contributions: MB, AM, CLM, and K-LL conceived the 
project and designed experimental details. MB performed hiPSC 
reprogramming together with CBJ. SC performed the single-cell 
electrophysiology measurements together with AW. LGJT did the 
MEA experiments. The targeting strategy was designed by RPD 
and MB, and experimental targeting work was done by MB and DW. 
DE provided the iViC^2. 5-hESC line. CD'A performed the FAGS 
sort and western blot experiments. JH did the HERG immunofluor- 
escence staining. The manuscript was written by MB, CLM, 
and AM. 

Conflict of interest 

Christine Mummery is co-founder and advisor of Pluriomics. The 
remaining authors declare that they have no conflict of interest. 



3172 The EMBO Journal VOL 32 | NO 24 | 2013 



©2013 European Molecular Biology Organization 



Isogenic pairs of LQT2 pluripotent stem cells 

M Bellin et al 



References 



Abyzov A, Mariani J, Palejev D, Zhang Y, Haney MS, Tomasini L, 
Ferrandino AF, Rosenberg Belmaker LA, Szekely A, Wilson M, 
Kocabas A, Calixto NE, Grigorenko EL, Huttner A, Chawarska K, 
Weissman S, Urban AE, Gerstein M, Vaccarino FM (2012) Somatic 
copy number mosaicism in human skin revealed by induced 
pluripotent stem cells. Nature 492: 438-442 

Asuri P, Bartel MA, Vazin T, Jang JH, Wong TB, Schaffer DV (2012) 
Directed evolution of adeno-associated virus for enhanced gene 
delivery and gene targeting in human pluripotent stem cells. Mol 
Ther 20: 329-338 

Axelrod FB, Liebes L, Gold-Von Simson G, Mendoza S, Mull J, 
Leyne M, Norcliffe-Kaufmann L, Kaufmann H, Slaugenhaupt SA 
(2011) Kinetin improves IKBKAP mRNA splicing in patients with 
familial dysautonomia. Pediatr Res 70: 480-483 

Barry PH, Lynch JW (1991) Liquid junction potentials and small cell 
effects in patch-clamp analysis. J Membr Biol 121: 101-117 

Bellin M, Marchetto MC, Gage FH, Mummery CL (2012) Induced 
pluripotent stem cells: the new patient? Nat Rev Mol Cell Biol 13: 
713-726 

Boulting GL, Kiskinis E, Croft GF, Amoroso MW, Oakley DH, 
Wainger BJ, Williams DJ, Kahler DJ, Yamaki M, Davidow L, 
Rodolfa CT, Dimos JT, Mikkilineni S, MacDermott AB, Woolf 
CJ, Henderson CE, Wichterle H, Eggan K (2011) A functionally 
characterized test set of human induced pluripotent stem cells. 
Nat Biotechnol 29: 279-286 

Braam SR, Tertoolen L, Casini S, Matsa E, Lu HR, Teisman A, 
Passier R, Denning C, Gallacher DJ, Towart R, Mummery CL 
(2013) Repolarization reserve determines drug responses in 
human pluripotent stem cell derived cardiomyocytes. Stem Cell 
Res 10: 48-56 

Chen X, Shen J, Prywes R (2002) The luminal domain of ATF6 
senses endoplasmic reticulum (ER) stress and causes transloca- 
tion of ATF6 from the ER to the Golgi. J Biol Chem 277: 
13045-13052 

Cherry AB, Daley GQ (2013) Reprogrammed cells for disease 
modeling and regenerative medicine. Annu Rev Med 64: 277-290 

Costa M, Dottori M, Sourris K, Jamshidi P, Hatzistavrou T, Davis R, 
Azzola L, Jackson S, Lim SM, Pera M, Elefanty AG, Stanley EG 
(2007) A method for genetic modification of human embryonic 
stem cells using electroporation. Nat Protoc 2: 792-796 

Crotti L, Celano G, Dagradi F, Schwartz PJ (2008) Congenital long 
QT syndrome. Orphanet J Rare Dis 3:18 

Davis RP, Casini S, van den Berg CW, Hoekstra M, Remme CA, 
Dambrot C, Salvatori D, Oostwaard DW, Wilde AA, Bezzina CR, 
Verkerk AO, Freund C, Mummery CL (2012) Cardiomyocytes 
derived from pluripotent stem cells recapitulate electrophysiolo- 
gical characteristics of an overlap syndrome of cardiac sodium 
channel disease. Circulation 125: 3079-3091 

Davis RP, Costa M, Grandela C, Holland AM, Hatzistavrou T, 
Micallef SJ, Li X, Goulburn AL, Azzola L, Elefanty AG, 
Stanley EG (2008) A protocol for removal of antibiotic resistance 
cassettes from human embryonic stem cells genetically modified 
by homologous recombination or transgenesis. Nat Protoc 3: 
1550-1558 

Davis RP, van den Berg CW, Casini S, Braam SR, Mummery CL 
(2011) Pluripotent stem cell models of cardiac disease and their 
implication for drug discovery and development. Trends Mol Med 
17: 475-484 

de Boer TP, van Veen TA, Jonsson MK, Kok BG, Metz CH, Sluijter JP, 
Doevendans PA, de Bakker JM, Goumans MJ, van der Heyden 
MA (2010) Human cardiomyocyte progenitor cell-derived cardio- 
myocytes display a maturated electrical phenotype. J Mol Cell 
Cardiol 48: 254-260 

Egashira T, Yuasa S, Suzuki T, Aizawa Y, Yamakawa H, Matsuhashi T, 
Ohno Y, Tohyama S, Okata S, Seki T, Kuroda Y, Yae K, 
Hashimoto H, Tanaka T, Hattori F, Sato T, Miyoshi S, Takatsuki S, 
Murata M, Kurokawa J et al (2012) Disease characterization using 
LQTS-specific induced pluripotent stem cells. Cardiovasc Res 95: 
419-429 

Ellgaard L, Helenius A (2003) Quality control in the endoplasmic 

reticulum. Nat Rev Mol Cell Biol 4: 181-191 
Elliott DA, Braam SR, Koutsis K, Ng ES, Jenny R, Lagerqvist EL, 

Biben C, Hatzistavrou T, Hirst CE, Yu QC, Skelton RJ, Ward-van 



Oostwaard D, Lim SM, Khammy 0, Li X, Hawes SM, Davis RP, 
Goulburn AL, Passier R, Prall OW et al (2011) NKX2-5(eGFP/w) 
hESCs for isolation of human cardiac progenitors and cardiomyo- 
cytes. Nat Methods 8: 1037-1040 

Fu J, Teucher M, Anastassiadis K, Skarnes W, Stewart AF (2010) A 
recombineering pipeline to make conditional targeting con- 
structs. Methods Enzymol 477: 125-144 

Giudicessi JR, Ackerman MJ (2013) Determinants of incomplete 
penetrance and variable expressivity in heritable cardiac arrhyth- 
mia syndromes. Transl Res 161: 1-14 

Gold-von Simson G, Goldberg JD, Rolnitzky LM, Mull J, Leyne M, 
Voustianiouk A, Slaugenhaupt SA, Axelrod FB (2009) Kinetin in 
familial dysautonomia carriers: implications for a new therapeu- 
tic strategy targeting mRNA splicing. Pediatr Res 65: 341-346 

Goldstein DB (2009) Common genetic variation and human traits. 
N Engl J Med 360: 1696-1698 

Gonzalez F, Boue S, Izpisua Belmonte JC (2011) Methods for making 
induced pluripotent stem cells: reprogramming a la carte. Nat Rev 
Genet 12: 231-242 

Gore A, Li Z, Fung HL, Young JE, Agarwal S, Antosiewicz-Bourget J, 
Canto I, Giorgetti A, Israel MA, Kiskinis E, Lee JH, Loh YH, 
Manos PD, Montserrat N, Panopoulos AD, Ruiz S, Wilbert ML, 
Yu J, Kirkness EF, Izpisua Belmonte JC et al. (2011) Somatic 
coding mutations in human induced pluripotent stem cells. 
Nature 471: 63-67 

Guo J, Zhang X, Hu Z, Zhuang Z, Zhu Z, Chen Z, Chen W, 
Zhao Z, Zhang C, Zhang Z (2012) A422T mutation in HERG 
potassium channel retained in ER is rescurable by pharmacologic 
or molecular chaperones. Biochem Biophys Res Commun 422: 
305-310 

Haze K, Yoshida H, Yanagi H, Yura T, Mori K (1999) Mammalian 
transcription factor ATF6 is synthesized as a transmembrane 
protein and activated by proteolysis in response to endoplasmic 
reticulum stress. Mol Biol Cell 10: 3787-3799 

Hetz C (2012) The unfolded protein response: controlling cell fate 
decisions under ER stress and beyond. Nat Rev Mol Cell Biol 13: 
89-102 

Hockemeyer D, Wang H, Kiani S, Lai CS, Gao Q, Cassady JP, 
Cost GJ, Zhang L, Santiago Y, Miller JC, Zeitler B, Cherone JM, 
Meng X, Hinkley SJ, Rebar EJ, Gregory PD, Urnov FD, Jaenisch R 

(2011) Genetic engineering of human pluripotent cells using TALE 
nucleases. Nat Biotechnol 29: 731-734 

Hoekstra M, Mummery CL, Wilde AA, Bezzina CR, Verkerk AO 

(2012) Induced pluripotent stem cell derived cardiomyocytes as 
models for cardiac arrhythmias. Front Physiol 3: 346 

Huffaker SJ, Chen J, Nicodemus KK, Sambataro F, Yang F, Mattay V, 
Lipska BK, Hyde TM, Song J, Rujescu D, Giegling I, Mayilyan K, 
Proust MJ, Soghoyan A, Caforio G, Callicott JH, Bertolino A, 
Meyer-Lindenberg A, Chang J, Ji Y et al (2009) A primate- 
specific, brain isoform of KCNH2 affects cortical physiology, 
cognition, neuronal repolarization and risk of schizophrenia. 
Nat Med 15: 509-518 

Itzhaki I, Maizels L, Huber I, Zwi-Dantsis L, Caspi 0, Winterstern A, 
Feldman 0, Gepstein A, Arbel G, Hammerman H, Boulos M, 
Gepstein L (2011) Modelling the long QT syndrome with induced 
pluripotent stem cells. Nature 471: 225-229 

Jung CB, Moretti A, Mederos y Schnitzler M, lop L, Storch U, Bellin M, 
Dorn T, Ruppenthal S, Pfeiffer S, Goedel A, Dirschinger RJ, 
Seyfarth M, Lam JT, Sinnecker D, Gudermann T, Lipp P, 
Laugwitz KL (2012) Dantrolene rescues arrhythmogenic RYR2 
defect in a patient-specific stem cell model of catecholaminergic 
polymorphic ventricular tachycardia. EMBO Mol Med 4: 180-191 

Khositseth A, Tester DJ, Will ML, Bell CM, Ackerman MJ (2004) 
Identification of a common genetic substrate underlying post- 
partum cardiac events in congenital long QT syndrome. Heart 
Rhythm 1: 60-64 

Kurien BT, Scofield RH (2006) Western blotting. Methods 38: 
283-293 

Lahti AL, Kujala VJ, Chapman H, Koivisto AP, Pekkanen-Mattila M, 
Kerkela E, Hyttinen J, Kontula K, Swan H, Conklin BR, Yamanaka S, 
Silvennoinen 0, Aalto-Setala K (2012) Model for long QT syndrome 
type 2 using human iPS cells demonstrates arrhythmogenic char- 
acteristics in cell culture. Dis Model Mech 5: 220-230 



©2013 European Molecular Biology Organization 



The EMBO Journal VOL 32 j NO 24 j 2013 3173 



Isogenic pairs of LQT2 pluripotent stem cells 

M Bellin et al 



Lee G, Papapetrou EP, Kim H, Chambers SM, Tomishima MJ, 
Fasano CA, Ganat YM, Menon J, Shimizu F, Viale A, Tabar V, 
Sadelain M, Studer L (2009) Modelling pathogenesis and treat- 
ment of familial dysautonomia using patient-specific iPSCs. 
Nature 461: 402-406 

Liang P, Lan F, Lee AS, Gong T, Sanchez-Freire V, Wang Y, Diecke S, 
Sallam K, Knowles JW, Wang PJ, Nguyen PK, Bers DM, Robbins RC, 
Wu JC (2013) Drug screening using a library of human induced 
pluripotent stem cell-derived cardiomyocytes reveals disease-spe- 
cific patterns of cardiotoxicity. Circulation 127: 1677-1691 

Ma J, Guo L, Fiene SJ, Anson BD, Thomson JA, Kamp TJ, Kolaja KL, 
Swanson BJ, January CT (2011) High purity human-induced 
pluripotent stem cell-derived cardiomyocytes: electrophysiologi- 
cal properties of action potentials and ionic currents. Am J Physiol 
Heart Circ Physiol 301: H2006-H2017 

Marchetto MC, Carromeu C, Acab A, Yu D, Yeo GW, Mu Y, Chen G, 
Gage FH, Muotri AR (2010) A model for neural development and 
treatment of Rett syndrome using human induced pluripotent 
stem cells. Cell 143: 527-539 

Matsa E, Dixon JE, Medway C, Georgiou 0, Patel MJ, Morgan K, 
Kemp PJ, Staniforth A, Mellor 1, Denning C (2013) Allele-specific 
RNA interference rescues the long-QT syndrome phenotype 
in human-induced pluripotency stem cell cardiomyocytes. 
Eur Heart J (advance online publication, 6 March 2013; 
doi:10.1093/eurheartj/eht067) 

Matsa E, Rajamohan D, Dick E, Young L, Mellor I, Staniforth A, 
Denning C (2011) Drug evaluation in cardiomyocytes derived 
from human induced pluripotent stem cells carrying a long QT 
syndrome type 2 mutation. Eur Heart J 32: 952-962 

McKernan KJ, Peckham HE, Costa GL, McLaughlin SF, Fu Y, 
Tsung EF, Clouser CR, Duncan C, Ichikawa JK, Lee CC, 
Zhang Z, Ranade SS, Dimalanta ET, Hyland FC, Sokolsky TD, 
Zhang L, Sheridan A, Fu H, Hendrickson CL, Li B et al (2009) 
Sequence and structural variation in a human genome uncovered 
by short-read, massively parallel ligation sequencing using 
two-base encoding. Genome Res 19: 1527-1541 

Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Each RP, 
Ingle C, Nisbett J, Guigo R, Dermitzakis ET (2010) 
Transcriptome genetics using second generation sequencing in a 
Caucasian population. Nature 464: 773-777 

Moretti A, Bellin M, Welling A, Jung CB, Lam JT, Bott-Flugel L, 
Dorn T, Goedel A, Hohnke C, Hofmann F, Seyfarth M, 
Sinnecker D, Schomig A, Laugwitz KL (2010) Patient-specific 
induced pluripotent stem-cell models for long-QT syndrome. 
N Engl J Med 363: 1397-1409 

Morita H, Wu J, Zipes DP (2008) The QT syndromes: long and short. 
Lancet 372: 750-763 

Moss AJ, Shimizu W, Wilde AA, Towbin JA, Zareba W, Robinson JL, 
Qi M, Vincent CM, Ackerman MJ, Kaufman ES, Hofman N, 
Seth R, Kamakura S, Miyamoto Y, Goldenberg I, Andrews ML, 
McNitt S (2007) Clinical aspects of type-1 long-QT syndrome by 
location, coding type, and biophysical function of mutations 
involving the KCNQl gene. Circulation 115: 2481-2489 

Moss AJ, Zareba W, Kaufman ES, Gartman E, Peterson DR, 
Benhorin J, Towbin JA, Keating MT, Priori SG, Schwartz PJ, 
Vincent CM, Robinson JL, Andrews ML, Feng C, Hall WJ, 
Medina A, Zhang L, Wang Z (2002) Increased risk of arrhythmic 
events in long-QT syndrome with mutations in the pore region of 
the human ether-a-go-go-related gene potassium channel. 
Circulation 105: 794-799 

Muller FJ, Schuldt BM, Williams R, Mason D, Altun G, 
Papapetrou EP, Danner S, Goldmann JE, Herbst A, Schmidt NO, 
Aldenhoff JB, Laurent LC, Loring JF (2011) A bioinformatic assay 
for pluripotency in human cells. Nat Methods 8: 315-317 

Mummery C, Ward-van Oostwaard D, Doevendans P, Spijker R, van 
den Brink S, Hassink R, van der Heyden M, Opthof T, Pera M, 
de la Riviere AB, Passier R, Tertoolen L (2003) Differentiation of 
human embryonic stem cells to cardiomyocytes: role of coculture 
with visceral endoderm-like cells. Circulation 107: 2733-2740 

Mummery CL, Zhang J, Ng ES, Elliott DA, Elefanty AG, Kamp TJ 
(2012) Differentiation of human embryonic stem cells and in- 
duced pluripotent stem cells to cardiomyocytes: a methods over- 
view. Circ Res 111: 344-358 

Narsinh KH, Plews J, Wu JC (2011) Comparison of human induced 
pluripotent and embryonic stem cells: fraternal or identical twins? 
Mol Ther 19: 635-638 



Nazor KL, Altun G, Lynch C, Tran H, Harness JV, Slavin I, 
Garitaonandia I, Muller FJ, Wang YC, Boscolo FS, Fakunle E, 
Dumevska B, Lee S, Park HS, Olee T, D'Lima DD, Semechkin R, 
Parast MM, Galat V, Laslett AL et al (2012) Recurrent variations in 
DNA methylation in human pluripotent stem cells and their 
differentiated derivatives. Cell Stem Cell 10: 620-634 

Panopoulos AD, Ruiz S, Izpisua Belmonte JC (2011) iPSCs: induced 
back to controversy. Cell Stem Cell 8: 347-348 

Pillozzi S, Brizzi MF, Balzi M, Crociani 0, Cherubini A, Guasti L, 
Bartolozzi B, Becchetti A, Wanke E, Bernabei PA, Olivotto M, 
Pegoraro L, Arcangeli A (2002) HERG potassium channels 
are constitutively expressed in primary human acute myeloid 
leukemias and regulate cell proliferation of normal and leukemic 
hemopoietic progenitors. Leukemia 16: 1791-1798 

Pini G, Scusa MF, Congiu L, Benincasa A, Morescalchi P, Bottiglioni I, 
Di Marco P, Borelli P, Bonuccelli U, Della-Chiesa A, Prina-Mello A, 
Tropea D (2012) IGFl as a potential treatment for Rett syndrome: 
safety assessment in six Rett patients. Autism Res Treat 2012: 679801 

Reinhardt P, Schmid B, Burbulla LF, Schondorf DC, Wagner L, 
Glatza M, Hoing S, Hargus G, Heck SA, Dhingra A, Wu G, 
Muller S, Brockmann K, Kluba T, Maisel M, Kruger R, Berg D, 
Tsytsyura Y, Thiel CS, Psathaki OE et al (2013) Genetic correction 
of a LRRK2 mutation in human iPSCs links Parkinsonian neuro- 
degeneration to ERK-dependent changes in gene expression. Cell 
Stem Cell 12: 354-367 

Rosati B, Marchetti P, Crociani 0, Lecchi M, Lupi R, Arcangeli A, 
Olivotto M, Wanke E (2000) Glucose- and arginine-induced 
insulin secretion by human pancreatic beta-cells: the role 
of HERG K( + ) channels in firing and release. FASEB J 14: 
2601-2610 

Rudy Y (2008) Molecular basis of cardiac action potential repolar- 
ization. Ann NY Acad Sci 1123: 113-118 

Sinnecker D, Goedel A, Laugwitz KL, Moretti A (2013) Induced 
pluripotent stem cell-derived cardiomyocytes: a versatile tool for 
arrhythmia research. Circ Res 112: 961-968 

Smith JL, McBride CM, Nataraj PS, Bartos DC, January CT, 
Dehsle BP (2011) Trafficking-deficient hERG K( + ) channels 
linked to long QT syndrome are regulated by a microtubule- 
dependent quality control compartment in the ER. Am J Physiol 
Cell Physiol 301: C75-C85 

Soldner F, Laganiere J, Cheng AW, Hockemeyer D, Gao Q, 
Alagappan R, Khurana V, Golbe LI, Myers RH, Lindquist S, 
Zhang L, Guschin D, Fong LK, Vu BJ, Meng X, Urnov FD, Rebar 
EJ, Gregory PD, Zhang HS, Jaenisch R (2011) Generation of 
isogenic pluripotent stem cells differing exclusively at two early 
onset Parkinson point mutations. Cell 146: 318-331 

Szuhai K, Tanke HJ (2006) COBRA: combined binary ratio labeling 
of nucleic-acid probes for multi-color fluorescence in situ hybri- 
dization karyotyping. Nat Protoc 1: 264-275 

Takahashi K, Okita K, Nakagawa M, Yamanaka S (2007a) Induction 
of pluripotent stem cells from fibroblast cultures. Nat Protoc 2: 
3081-3089 

Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, 
Yamanaka S (2007b) Induction of pluripotent stem cells from 
adult human fibroblasts by defined factors. Cell 131: 861-872 

Tan HL, Bardai A, Shimizu W, Moss AJ, Schulze-Bahr E, Noda T, 
Wilde AA (2006) Genotype-specific onset of arrhythmias in 
congenital long-QT syndrome: possible therapy implications. 
Circulation 114: 2096-2103 

Tester DJ, Will ML, Haglund CM, Ackerman MJ (2005) 
Compendium of cardiac channel mutations in 541 consecutive 
unrelated patients referred for long QT syndrome genetic testing. 
Heart Rhythm 2: 507-517 

Vandenberg JI, Perry MD, Perrin MJ, Mann SA, Ke Y, Hill AP (2012) 
hERG K( + ) channels: structure, function, and clinical signifi- 
cance. Physiol Rev 92: 1393-1478 

Wang Y, Huang X, Zhou J, Yang X, Li D, Mao H, Sun HH, Liu N, 
Lian J (2012) Trafficking-deficient G572R-hERG and E637K-hERG 
activate stress and clearance pathways in endoplasmic reticulum. 
PLoS One 7: e29885 

Yazawa M, Hsueh B, Jia X, Pasca AM, Bernstein JA, Hallmayer J, 
Dolmetsch RE (2011) Using induced pluripotent stem cells to 
investigate cardiac phenotypes in Timothy syndrome. Nature 471 : 
230-234 

Yu J, Vodyanik MA, Smuga-Otto K, Antosiewicz-Bourget J, Frane JL, 
Tian S, Nie J, Jonsdottir GA, Ruotti V, Stewart R, Slukvin II, 



3174 The EMBO Journal VOL 32 | NO 24 | 2013 



©2013 European Molecular Biology Organization 



Isogenic pairs of LQT2 pluripotent stem cells 

M Bellin et al 



Thomson JA (2007) Induced pluripotent stem cell lines derived 
from human somatic cells. Science 318: 1917-1920 
Zou J, Maeder ML, Mali P, Pruett-Miller SM, Thibodeau-Beganny S, 
Chou BK, Chen G, Ye Z, Park IH, Daley GQ, Porteus MH, Joung 
JK, Cheng L (2009) Gene targeting of a disease-related gene in 
human induced pluripotent stem and embryonic stem cells. Cell 
Stem Cell 5: 97-110 



[Z^ ® I The EMBO Journal is published by Nature 
^^i^^ms^m Publishing Group on behalf of the European 
Molecular Biology Organization. This article is licensed 
under a Creative Commons Attribution 3.0 Unported 
Licence. To view a copy of this licence visit http://creative 
commons.org/licenses/by/3 .0/. 



©2013 European Molecular Biology Organization 



The EMBO Journal VOL 32 j NO 24 j 2013 3175 



